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A fermion dark matter candidate with a relic abundance set by annihilation through a pseu¬ 
doscalar can evade constraints from direct detection experiments. We present simplified models 
that realize this fact by coupling a fermion dark sector to a two-Higgs doublet model. These mod¬ 
els are generalizations of mixed bino-Higgsino dark matter in the MSSM, with more freedom in 
the couplings and scalar spectra. Annihilation near a pseudoscalar resonance allows a significant 
amount of parameter space for thermal relic dark matter compared to singlet-doublet dark matter, 
in which the fermions couple only to the SM Higgs doublet. In a general two-Higgs doublet model, 
there is also freedom for the pseudoscalar to be relatively light and it is possible to obtain thermal 
relic dark matter candidates even below 100 GeV. In particular, we find ample room to obtain dark 
matter with mass around 50 GeV and fitting the Galactic Center excess in gamma-rays. This region 
of parameter space can be probed by LHC searches for heavy pseudoscalars or electroweakinos, and 
possibly by other new collider signals. 


I. INTRODUCTION 

Is weak-scale dark matter (DM) still a viable sce¬ 
nario? At face value, a DM candidate of mass 10- 
1000 GeV with weak-scale interactions can easily have 
a thermal relic abundance matching the observed value 
of Ulcdmh'^ ~ 0.1199 ± 0.0027 [1]. However, the annihi¬ 
lation process is generically related by crossing symme¬ 
try to interactions in direct detection experiments and 
at colliders, both of which are becoming increasingly re¬ 
strictive for this range of DM masses. To satisfy these 
constraints requires more ingredients or tunings in many 
existing models. 

One simple and attractive possibility is that the DM 
interaction with Standard Model (SM) particles is sup¬ 
pressed in the non-relativistic limit in the t-channel (di¬ 
rect detection) while preserving a weak-scale cross section 
in the s-channel (for a thermal relic). This requirement is 
satisfied by a pseudoscalar mediator coupling to fermion 
DM and to SM fermions: 

^ D y^Axii^x + . (1) 

Because the pseudoscalar interaction breaks chiral sym¬ 
metry, the coupling of A with SM fermions is generically 
proportional to the SM Yukawa couplings yf (A/ oc yf). 
Therefore, since the pseudoscalar has larger couplings to 
third-generation quarks and leptons, collider constraints 
are typically weaker. Integrating out the pseudoscalar 
gives the dimension six operators 

^(A7^x){/7V) • (2) 

rriA 

These contact operators have been considered in Refs. [2- 
7] , which motivated collider signals with MET and a sin¬ 
gle jet or 6-jet as a new search channel for DM. 

In order to provide a concrete but simple realization 
of this interaction, the approach taken in this paper is to 


build a simplified model of DM coupled to a new pseu¬ 
doscalar mediator. The philosophy is to add a minimal 
set of new matter fields with renormalizable and gauge- 
invariant couplings [8]. Our models also provide a UV- 
completion (at least at scales relevant to the LHG) of the 
type of contact interactions discussed above. 

Since the SM does not contain a fundamental pseu¬ 
doscalar, we focus our study on two-Higgs doublet models 
(2HDM). We note that the simplest cases of Higgs-portal 
models connecting DM to the SM through Higgs media¬ 
tion are highly constrained [9-11]. In particular, existing 
models have studied singlet-doublet fermion sectors cou¬ 
pled to the SM Higgs [12-14], while our work generalizes 
this to 2HDMs. Within the 2HDM framework, the pseu¬ 
doscalar interaction with DM carries unavoidable inter¬ 
actions of DM to the CP-even scalars as well. However, 
compared to the usual Higgs portal models, in a general 
2HDM, it is possible to obtain parametrically larger cou¬ 
plings of DM to the pseudoscalar (or heavy Higgs) than 
to the light Higgs, in this way alleviating the various con¬ 
straints. 

The possibility of resonant annihilation through the 
pseudoscalar also permits a more sizable parameter space 
for DM to be a thermal relic. As we will show, much of 
the open parameter space is in the region where ^ 
mAj^. This is the analog of the so-called “A-funnel” 
region in the Minimal Supersymmetric Standard Model 
(MSSM). Still, it is important to note that, even if the 
DM mass is not tuned to be very close to m^/2, we can 
obtain a thermal DM and satisfy the constraints from 
direct detection and from the LHC. 

Another motivation to consider models with DM 
mainly annihilating through the A-funnel comes from the 
intriguing results for DM annihilation in the gamma-ray 
sky [15-21]. Studies of the Fermi gamma-ray data show 
an excess diffuse component around the Galactic Cen¬ 
ter (GC). Although these diffuse gamma-rays may have 



an astrophysical origin, it is interesting to note that their 
spectrum and morphology is consistent with expectations 
for gamma-rays from DM annihilation. The spectrum 
peaks at around 2-3 GeV, which suggests a DM candi¬ 
date with mass smaller than around 100 GeV, and the 
signal follows an NFW profile. Many theoretical stud¬ 
ies have explored simplihed scenarios [22-26], as well as 
more UV-complete models, for the Galactic Center Ex¬ 
cess (GCE). In particular, light DM with mass around 
~ 50 GeV annihilating mainly to bottom quarks and tau 
leptons has identified as giving a good fit to the data. 
However, a more recent Fermi analysis of the GC re¬ 
gion [27] indicates the spectrum of the diffuse component 
is harder than previously thought, allowing for reasonable 
fits with heavier DM candidates annihilating to di-boson 
final states [28-30]. 

A number of previous works have identified the 
pseudoscalar case as a promising candidate to ht the 
GCE [31-41]. In a number of these models, in order 
to satisfy various experimental constraints, a new light 
pseudoscalar which mixes with the pseudoscalar in the 
2HDM is introduced, while the new scalars in the 2HDM 
itself are relatively decoupled [32-34, 42]. The new fea¬ 
ture of our model is that we rely purely on the pseu¬ 
doscalar of the 2HDM, but consider more general 2HDMs 
and dark sector spectra. 

The paper is organized as follows. In Section II we 
review aspects of general 2HDMs and we fix our nota¬ 
tion. Furthermore, we motivate the scalar spectra con¬ 
sidered in the following sections. Section III discusses the 
model-independent collider constraints on the simplihed 
model, covering heavy Higgs searches and searches with 
MET. We then turn to specihc models, beginning with a 
brief discussion on extended Higgs portal models for both 
scalar and fermion DM in Section IV. Our main analysis 
is contained in Section V for fermion DM coupled to a 
2HDM. We consider a range of constraints on the model, 
and present viable parameter space both for heavier DM 
(^150 GeV) in Section V G and specihcally for the GCE 
in Section V D. We reserve Section VI for our conclusions. 
In the appendices, we elaborate on 2HDM benchmarks 
for the GCE (Appendix A) and give analytic formulae 
for Higgs couplings to DM in our model (Appendix B). 


II. TWO HIGGS DOUBLETS 


The most general renormalizable Higgs potential of a 
2HDM can be written as 


V = 


( 3 ) 


+A3($i$,)($t $„) + A4($^4>„)(4>i$d) 
+A 6 ($d$^)$d^d + A 7 ($i$n)$t$„ + h.c.l , 


with and Higgs doublets with hypercharge -I-I/2. 
The Higgs helds can be parametrized as 


= 


= 
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Assuming CP conservation, the mass eigenstates are 
given by 



/cos a — sina\ / hu\ 
l^sina cos a J \hdj ’ 

/ sin/3 cos/3 A /a^A 
[cosP -s\n(3j yadJ ’ 

/sin/3 cos/ 3A/H±A 
l^cos /3 — sin /3y \HiJ ’ 


( 5 ) 


where we dehne the basis-dependent ratio 
tan/3 = Vu/vd ^ and = 246 GeV. 

and are Goldstone bosons, A the pseudoscalar, 
the charged Higgs, and h, H the light and heavy CP-even 
scalar, respectively. In the following, we will choose to 
work in the alignment limit, sin(/3 — a) = I, where the 
h couplings are SM-like. 

In the MSSM, the masses of the heavy scalars are clus¬ 
tered around a similar scale, with splittings arising only 
from small D-terms. However, in a general 2HDM frame¬ 
work, we have more freedom to get more sizable split¬ 
tings. As a result, we can split the pseudoscalar mass 
such that mA < mn ~ 'mH+, as needed for a model that 
can ht the GCE (see Sec. VD). In particular, we can 
write the mass eigenvalues as functions of the quartic 
couplings in Eq. (3): 


= ^(As - A 4 ) ( 6 ) 

and for uia 3^ Kv and tan/3 ^ I, 


— m\ ~ + A 2 ~ 2 (A 3 -I- A 4 — A 5 )), (7) 

mi ~ ^(Ai -I- A 2 + 2(A3 -I- A 4 + A 5 -I- 2 A 67 )), ( 8 ) 


while for mA 3^ AiU and tan/3 S> I, 

mi — m\ ~ Asu^, m\ ~ . (9) 


where A 67 = Ag -I- A 7 . From these expressions, it is clear 
that a splitting as large as ^100 GeV between the pseu¬ 
doscalar and the charged Higgs and between the pseu¬ 
doscalar and the heavy Higgs H can be obtained for 
A4,A5 = 0(I). 

Electroweak precision measurements, and in particu¬ 
lar the p parameter, can give important constraints on 
spectra with large splittings. A notable exception are 


^ In the following, we will restrict our attention to Type-II 2HDMs, 
for which tan /3 is a well defined basis-independent quantity [43] . 
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alignment models with niH ~ , but with an arbi¬ 

trary large mass splitting between the pseudoscalar and 
the charged Higgs, which leads to a very small correction 
to the p parameter. Further constraints arise from the re¬ 
quirement of vacuum stability and perturbativity of the 
quartic couplings. In Appendix A, we will give examples 
of viable quartic couplings that are able to produce the 
spectra that we use in Sec. VD, including the right mass 
for the SM-like Higgs boson discovered at the LHC. 

Finally, additional constraints on the scalar spectrum 
come from flavor transitions. In general, both Higgs dou¬ 
blets can couple to up and down quarks, as well as to 
leptons, with generic Yukawa couplings. However, the 
most general Yukawa couplings lead to excessively large 
contributions to flavor changing neutral transitions. This 
has led to the various well-known “Types” of 2HDMs, in 
which different discrete Z 2 symmetries are imposed; the 
two Higgs doublets have different Z 2 charge, thus for¬ 
bidding some of the Yukawa couplings. This condition 
also reduces the number of free parameters in Eq. (3), 
since the Z 2 symmetry demands Ag = Ay = 0^. In the 
following, we will focus on a Type-H 2HDM, for which 
the Z 2 symmetry allows the doublet to only couple to 
right-handed up quarks and the doublet <I>d to only cou¬ 
ple to right-handed down quarks and leptons. This type 
of 2HDM is interesting, since the MSSM at tree level is 
a particular limit of Type-H 2HDMs. 


III. COLLIDER SIGNALS 


In this section we present model-independent collider 
bounds on the new Higgs sector coupled to DM. We will 
keep the discussion as general as possible, considering 
a single new pseudoscalar and dark matter particle with 
coupling as in Eq. (1). Collider signals of simplified mod¬ 
els of dark matter coupled to new scalars have also been 
discussed recently in [44, 45]. For the set of models we 
consider in this paper, there will also be a rich set of 
collider signals associated with new charged states, such 
as the charged Higgs and the additional charged fermion 
states. We will discuss these in Section V. 


A. Present bounds 


In a Type-H 2HDM, the A coupling to down-type 
quarks, and charged leptons, is given by y/tan/3, where 
yf = rrif/v^ with ruf the mass of the corresponding 
fermion. Assuming a Majorana fermion DM x with in¬ 
teraction as in Eq. (I), the invisible decay width of the 


^ We work under the assumption that the Z 2 symmetry is softly 
broken by the mass term + h.c.. 



FIG. 1. Comparison of present LHC bounds for a heavy pseu¬ 
doscalar in the visible and invisible decay channels, in each 
case assuming the maximum possible branching ratio arising 
in a Type-II 2HDM (solid lines). Upper limits on tan/J are 
shown as a function of the new scalar mass ruA- Note that for 
A ^ bb only 7 TeV limits are available, while all the others 
searches have been performed with the full 8 TeV, 20/fb data 
set. We also show extrapolated limits for a 14 TeV LHC with 
300/fb (dashed lines), as explained in Section HIB. 


heavy pseudoscalar is given by 


r(A ^xx) = y 


2 rriA 
X 4,^ 



while the visible decay width is 


r(A ^ //) 


ricyj tan^ j3 


rriA 

Stt 



( 10 ) 


( 11 ) 


where ric is the number of colors of the final state SM 
fermion. Depending on the values of y^ and tan/?, dif¬ 
ferent LHC search channels can play a role. 

In Fig. I we compare the constraints on the visible 
and invisible A decay from current LHC searches. In the 
same figure, we also give 14 TeV projections for each of 
these searches, as discussed in Section HIB. In each case, 
we assume the maximal possible branching ratio in order 
to show the optimal sensitivity that could be achieved 
in each channel. More precisely, in the case of the pseu¬ 
doscalar decaying to bottom quarks and r leptons, we 
assume yx = ^ and in the case of the pseudoscalar decay¬ 
ing invisibly, we assume BR(A —>■ xx) 1- For simplic¬ 
ity, results are shown in the narrow width limit for the 
pseudoscalar. 

In the absence of invisible decays, in a Type-H 2HDM 
at large tan /?, the branching ratio of A to rr is about 
10% and the one to bb about 90%. The most stringent 
LHC constraints on A come from heavy MSSM Higgs 
boson searches. The main production modes of A are 
in association with 5-quarks and in gluon fusion with a 
heavy quark loop. For A —)■ bb, we use the 7 TeV anal¬ 
ysis from CMS [46]. For A —>■ tt, we take the 8 TeV 
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FIG. 2. Present and projected collider reach on the pseu¬ 
doscalar couplings for a specific choice of DM and pseu¬ 
doscalar masses. The parameters chosen are those favored 
for explaining the GCE, as discussed in Section VD. We also 
show the 8 TeV limits from the ti -|- fix search [53] . 


results from ATLAS [47], which gives slightly stronger 
limits at high niA compared to the CMS analysis [48] . In 
order to derive conservative constraints, we assume that 
the heavy CP-even Higgs is somewhat heavier than the 
pseudoscalar, so that only the pseudoscalar contributes 
to the invisible or visible signature. This is different than 
the typical models used by the LHC experimental col¬ 
laborations to interpret heavy Higgs searches [46-48], in 
which mA = rnn and both the pseudoscalar and heavy 
CP-even Higgs contribute to the signature. 

The XX decay mode gives rise to a b+^T (mono-6) 
signal in 5-associated production. To derive constraints, 
we simulate the signal with MadGraphS [49], PYTHIA 
[50] and DELPHES [51], and compare with the signal 
region SRI in the ATLAS analysis [52]. As Fig. 1 shows, 
the TT channel gives the strongest constraint. This is 
the only channel for which a reconstruction of ttia with 
ttitt is possible, and furthermore can have smaller back¬ 
grounds if a leptonic r decay is considered. 

Fig. 1 assumed maximal branching ratios in each 
case. The relative importance of the constraints is mod¬ 
ified when both invisible and visible decays have sizable 
branching ratios. Next, we consider the interplay of these 
searches, as well as the requirement of a thermal relic. 
For concreteness, we focus on DM with mass = 50 
GeV, where the dominant annihilation channel is through 
a pseudoscalar of mass ttia = 160 GeV, such as that fa¬ 
vored by the GCE (see Section VD). 

The constraints on the two free parameters and 
tan/3 are shown in Fig. 2. For both visible and invisi¬ 
ble decays in 6-associated production, a minimum tan/3 
can be reached since the production mechanism relies on 
the Abb coupling, which is tan/3 enhanced. Furthermore, 
as expected, the mono-6 search covers the large region. 


which is difficult to probe with the 666 search. In Fig. 2, 
we also show limits at low tan /3 coming from the search 
for a pseudoscalar decaying invisibly and produced in as¬ 
sociation with top quarks [53]. This search is able to 
exclude values of tan/3 ~ 0.5, above which the produc¬ 
tion cross section becomes too small^. From the figure, 
we can see that the tt + and 6tt signatures are com¬ 
plementary to probe the region of parameter space where 
DM is a thermal relic. The b + l^x search, instead, is not 
effective in this region. 

Finally, in this context, it is interesting to understand 
where the SM-like Higgs boson discovered by the LHC is 
placed. The most stringent direct bound on the branch¬ 
ing ratio of the Higgs decaying invisibly comes from the 
CMS analysis [54] that combines the searches for a Higgs 
produced in vector-boson-fusion and for a Higgs pro¬ 
duced in association with a Z boson. This analysis is 
able to put a bound on the branching ratio into invisible 
at 58% at the 2 (t level. Additional bounds come from 
searches for a Higgs produced in association with tops 
and decaying invisibly. Following the results in Ref. [53], 
the branching ratio into invisible is bounded at 40%. Us¬ 
ing this latter bound, the SM Higgs coupling to DM is 
then constrained to be about 0.02, with which a thermal 
relic is possible in only-Higgs mediated models when 
is very close to 771^/2 [55]. This is a much stronger bound 
than we get from the h + ^x signature, which is only able 
to constrain at best a coupling of ~ 1 for the same 
scalar mass. 


B. LHC Run 2 projections 

In Figs. 1-2, we also show projections for the 14 TeV 
LHC with 300/fb integrated luminosity. 

For the mono-6 limits, we simulate events with the cuts 
in [2], taking advantage of the background simulation al¬ 
ready done there. A full projection for the visible chan¬ 
nels requires modeling of the gluon fusion to A processes, 
as well as the SM 666, tt and 6rr backgrounds. A de¬ 
tailed collider study is beyond the scope of the current 
paper. We present a simple estimate for these channels. 
We extrapolate both signal and background based on pdf 
luminosity, assuming for the signal on-shell production 
of a new heavy mediator. Our method is similar to that 
in [56], but we have also to include the 6—quark PDF in 
our analysis. Combining these results gives an estimate 
for the improvement in S/Vb, which we use to deter¬ 
mine the improvement in the tan /3 bound for each given 
niA- For this reason, in Fig. 1 the projection for the 666 
search is only shown up to itla = 350 GeV. In all cases. 


® In principle, invisible h decay should also be included which 
would strengthen constraints, but the relation of its invisible de¬ 
cay to the A invisible decay is more model-dependent, so we do 
not include it in this figure. 
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for A production, we consider only 6-associated produc¬ 
tion. The gluon fusion channel will become important at 
low values of tan /?. For the values of tan j3 probed at 8 
TeV, we have checked that gluon fusion is subdominant 
if compared to 6-associated production, with the explicit 
cross section limits on gg —>■ bbA given in Refs. [47, 48]. 
Therefore, the projections for 14 TeV can be regarded 
as conservative, in that they ignore the additional con¬ 
straining power from gluon fusion. 

Two comments on the figures are in order: projections 
of the bbb and brr channels show that, as expected, the 
bound on tan /3 from the A ^ bb decay will be weaker 
than the bound from A ^ tt. Still, the difference in the 
reach of the two channels is not that large, especially at 
low values of m^, motivating the search for A ^ bb with 
additional data beyond the 5/fb 7 TeV set. Additionally, 
as we can see from Fig. 2, improved constraints on the 
bTT (bbb) and tt+ MET channels from LHC run 2 could 
help cover much of the region of parameter space where 
DM is a thermal relic. 


IV. EXTENDED HIGGS PORTAL MODELS 


Models of scalar DM coupled to an extended Higgs 
sector have been studied extensively in the literature, in¬ 
cluding the case with a 2HDM Higgs content [57-61]. 
Especially for light DM, this is a much more constrained 
scenario than for fermion DM, as we review in the fol¬ 
lowing. 

If DM is assumed to be a real scalar, S, and a gauge 
singlet under the SM, the lowest dimension gauge invari¬ 
ant operator that allows direct couplings to the SM Higgs 
is [62] 

/:dm D ^ (^2) 


where d* is the SM Higgs doublet. However, for DM mass 
ms ^ mhl2 the large values of Ag needed to sufficiently 
deplete the abundance of S before freeze-out are in strong 
conflict with constraints from the invisible width of the 
SM Higgs. For heavier masses, aside from a small win¬ 
dow of viable parameter space very close to the Higgs 
resonance, LUX [63] rules out thermal DM for ms < 100 
GeV [55] (see also [64], for a recent analysis). However, 
this does not have to be the case if the model is slightly 
extended beyond the simple Lagrangian of Eq. (12). 

A concrete example is the singlet scalar extension of 
a 2HDM. Once again, the singlet scalar S is the DM 
candidate, but it now possesses interactions with ^d,u 
through the terms 


£dm a 


Ad^’d^'d + -k Ad„(4’]i4>„ -k h.c.) 


The introduction of a second Higgs doublet allows S to 
annihilate through an s-channel H as well. For CP- 
conserving interactions, annihilation through A is not al¬ 
lowed. As a result, it is possible to slightly uncouple the 
annihilation rate from the coupling of the SM Higgs with 


DM and thus ease the tension with the invisible Higgs 
and LUX constraints to some degree. However, doing 
so can require tunings in the couplings above, in order 
to sufficiently suppress the interaction strength with the 
SM Higgs. Furthermore, since the heavy H also leads to 
spin-independent (SI) nucleon scattering, thermal DM 
below 100 GeV is still very restricted by direct detection 
constraints. More detailed analysis of the model can be 
found in Ref. [57]. 

Another option is to consider fermion DM. We begin by 
summarizing the constraints on a model with the Higgs 
sector given just by the SM Higgs. Light thermal DM is 
particularly difficult to achieve. In particular, a strictly 
weakly interacting fermion DM candidate is in strong 
tension with the null results of current direct detection 
experiments. For example, a DM candidate like a 4th 
generation heavy Dirac neutrino has Z couplings that 
would lead to a nucleon elastic scattering rate many or¬ 
ders of magnitude beyond what is currently allowed by 
LUX (see e.g. [55]). 

A fix for this issue is to introduce a new gauge singlet 
fermion that is allowed to mix with the active compo¬ 
nents of the DM after electroweak symmetry breaking. 
This can be explicitly realized by coupling the sterile 
and active components through renormalizable interac¬ 
tions involving the SM Higgs. Dubbed “singlet-doublet” 
DM [12-14], the simplest possibility is to introduce a 
Majorana gauge singlet and a vector-like pair of SU(2) 
doublets. Since the mixing of the DM candidate orig¬ 
inates from Yukawa interactions with the SM Higgs, 
such scenarios generically are quite constrained by spin- 
independent direct detection limits, and for light DM 
this usually necessitates living somewhat near “blind- 
spots” of the theory, where the coupling to the SM Higgs 
is strongly suppressed. Furthermore, since annihilation 
through the SM Higgs is velocity suppressed, to accom¬ 
modate a thermal relic, DM annihilations involving gauge 
bosons need to be the dominant channel, and for DM 
lighter than ~ 80 GeV, this demands a strong coupling 
to the Z. However, precision electroweak measurements 
near the .Z-pole strongly constrain such a coupling when 
m^ < mzl2 [65]. Hence, additional annihilation media¬ 
tors would greatly benefit the prospects for light thermal 
DM. 

In light of this reasoning, in the next section we will 
study a model that adds a second Higgs doublet to this 
simple fermion DM scenario. The singlet-doublet case 
discussed above is one particular limit of this model. As 
we will see, the extension to a 2HDM vastly alleviates 
the tension of < 0(100) GeV singlet-doublet DM with 
current experimental constraints. 


V. FERMION DM IN A 2HDM 
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We define a simple UV-complete extension of the SM, 
whose scalar sector is that of a Type-H 2HDM, and whose 
DM sector consists of one gauge singlet Weyl fermion. S', 




Field 

Charges 

Spin 


(1,2, +1/2) 

0 


(1,2, +1/2) 

0 

S 

(1,1,0) 

1/2 

Di 

(1,2, -1/2) 

1/2 

D2 

(1,2, +1/2) 

1/2 


TABLE I. Field content of our model with fermion DM. The 
given charges are for the SM gauge groups SU{S) x SU{2) x 


Name 

Model 

du 

4>1 = '1>2 = <l>u 

ud 

= 4>u, 4>2 = 

dd 

■I?! = "Fd, $2 = ^d 

UU 

4>1 = 4>2 = 4’u 


TABLE 11. Naming scheme for Eq. (13). 


since we are freed from restrictions such as the Yukawa 
couplings being fixed by the gauge sector, and compared 
to singlino-Higgsino DM in the NMSSM, the model re¬ 
quires fewer new degrees of freedom^. For simplicity, we 
take all couplings to be real, and work in the convention 
Ms, Mjj > 0. Keeping the sign of Mu and Ms fixed, we 
see that we can parity transform just S —> —S or both 
Di —)■ —Di and D 2 —t — D 2 . Either of these choices in 
field re-definitions results in simultaneously flipping the 
sign of j/i and 2/2 ■ However, their relative sign remains 
unchanged. Therefore, the only physical sign in our cou¬ 
plings is that of y 2 /yi- Hence, we will often trade the 
couplings 2/1 , 2/2 in favor of y and tan 6 defined by 

2/1 = 2 /cosd , 2/2 = ysind . (14) 

The DM doublets are parametrized as 



and two Weyl SU{2) doublets Di and D 2 , oppositely 
charged under C/(l)y, as shown in Table I. We assume 
that the fermionic fields of the dark sector are Z 2 odd, 
ensuring the stability of the lightest state. 

The DM sector Lagrangian has the terms 

>Cdm D — -MsS^ — MdDiD2 

-2/i5Di$i-2/25$^D2+h.c. , (13) 


where ^ 1 ^ 2 , E’ 1,2 are the neutral, charged components of 
the doublets, respectively. Ei and E 2 combine to form 
an electrically charged Dirac fermion of mass Md. To 
avoid constraints coming from chargino searches at LEP, 
we only consider values of Md ^ 100 GeV [66]. 

After electroweak symmetry breaking, one ends up 
with three Majorana mass eigenstates, in general all mix¬ 
tures of the singlet S and the neutral components 
and V 2 ) of the SU{2) doublets. In the (S', 1 ^ 1 , 2 ^ 2 ) basis, 
the neutral mass matrix for the fermionic dark sector is 


where 2-component Weyl and SU(2) indices are implied. 
The notation $ 1,2 indicates that we will allow different 
permutations for which Higgs doublets (down or up-type) 
couple to Di and D 2 ■ From here forward, we will use the 
naming scheme where $1 = $£;, $2 = will be called 
the “du” model, and similarly for other choices of $1 
and $ 2 . The naming schemes are given in Table 11. Note 
that in Eq. (13), we have made the simplifying assump¬ 
tion that £>1 and D 2 each only couple to a single Higgs 
doublet. In the case of the dd and uu models, this can 
be enforced by the Z 2 symmetry of the 2HDM scalar sec¬ 
tor (see Sec. H), that is assumed to be broken only by 
mass terms. This is not the case for the du and ud mod¬ 
els. Still, our qualitative conclusions would not change 
significantly if all possible Yukawas are allowed to be of 
comparable strength. 

The situation where only the SM Higgs doublet ap¬ 
pears in the interactions of the Lagrangian (13) is the 
singlet-doublet model discussed in Section IV and can 
be identified as the limiting case of tan ,5 1 in the uu 

model. In this paper, we do not perform an analysis of 
this model, as it has already been covered in great detail 
in [13] and [14]. Furthermore, since the du and ud model 
can be mapped into each other by a simple switch of 
2/1 ^ 2 / 2 , only one systematic study of the two is needed. 

Although the du model embodies the spirit of bino- 
Higgsino DM in the MSSM, it allows much more freedom 


Afneutral 


^ Ms ^yivi j^y2V2'^ 
■^yivi 0 Md 

\^:^y2V2 Md 0 j 


(16) 


where vi ^2 are the VEVs of 4>i_2. The lightest mass eigen¬ 
state of Mneutrai will be the stable Majorana DM candi¬ 
date, which we write as a 2-component Weyl fermion Xa- 
We will denote the Xa composition in terms of the gauge 
eigenstates as 


Xa = + N^V2c. ■ (17) 


In most of the viable parameter space for thermal DM, 
especially for light DM < (!)(100) GeV, it is usually 


^ Bino-Higgsino mixing of the MSSM can be related through the 
identifications y -H- g', tan0 -H- —1, Ms ^ Mi, and Mf) -H- /i, 
where g' is the hypercharge gauge coupling, Mi is the bino soft 
SUSY breaking mass, and y is the Higgsino mass term. For 
singlino-Higgsino mixing of the NMSSM, the appropriate iden¬ 
tifications are y -H- ^/2X, tanO -H- -^1, Ms ^ 2acs, and 

Mj:) •<-> As, where A is the Yukawa coupling for the trilinear 
singlet-Higgs interaction of the superpotential, y.' is the super- 
symmetric mass term for the singlet, k is the Yukawa trilinear 
self-interaction for the singlet, and s is the VEV of the singlet 
scalar. 
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the case that the x is singlet-like and ^ Ms Md- 
This is because a large doublet component is strongly 
constrained by direct detection. Therefore, it is useful 
to write down the approximate form of Ng 2 when both 
Ms and are much smaller than Mo. In this limit, 

one finds 


iV| « 1 , Nf 


— /yx ~_L 

72 Md ’ 2 V2Md ' 


(18) 


A. Interactions 


We present the relevant DM interactions in this sec¬ 
tion. For the rest of the analysis, we will give our ana¬ 
lytic expressions in terms of 4-component notation, since 
they can be written in a more compact form. We write 
the charged state as E, a Dirac fermion of electric charge 
Q = —1 and mass Mo- 

The interactions of x with h, H, A, and are 

C D XxhhxX + XxhHxx + XxaAxii^x 
+ \H^X E + h.c.]. (19) 

The full forms for these couplings are given in Ap¬ 
pendix B. Note that in their explicit analytic form in 
Eqs. (B6)-(B9), the Higgs couplings are proportional to 
reiterating that the only physical sign is that of tan 9. 
All of these couplings rely on the singlet-doublet mix¬ 
ing. Also from Eqs. (B6)-(B9), it is important to note 
that, as long as Di and D 2 couple to the same Higgs 
doublet (i.e. the dd and uu models), there is a generic 
“blind-spot” of the theory that is independent of tan/I: 
for {m^ + Mo sin 2d) ~ 0 the couplings to the CP-even 
scalars identically vanish. We will often be considering 
the case where rrix Mo, since a large doublet com¬ 
ponent is strongly constrained by direct detection. In 
this limit, the blind-spot in these scenarios will necessi¬ 
tate |tand| <C 1 or |tand| ^ 1, along with tand < 0. 
Furthermore, for rux Mo, it is particularly simple to 
write the coupling of x to the neutral Higgs bosons as 


Xxh 

«ilF«^2±/v2%) 


\h 





(20) 


where are the projections of the physical Higgses 

into the gauge eigenstates, as defined in Eq. (Bl), and 
are functions of tan/3. 

Through mixing with vi and V 2 , X can also have inter¬ 
actions with the electroweak gauge bosons of the form 

C D gxzZ^xi^x^x + XT'" {g+v + 9+al^) E + h.c.] 

( 21 ) 

such that 


9xz = 


4r 


(/Vf)" - 


8c,„Mfy 


r 2 2 2 

[ViVi - 92^: 


9+v — 
9+a = 


272 

-9 

272 


(N^ + N}) 


(iVf - N}) 


^{yivi + y2V2) 

IMo ~ ’ 


( 22 ) 


where the latter approximate expressions hold for ^ 
Ms Mo- g is the 51/(2) gauge coupling and Cyj is 
the cosine of the Weinberg angle. It is apparent that, 
for Ms Mo, gxZ is small, and, furthermore, even for 
large 51/(2) doublet mixing, g^z can vanish identically 
if Nf = , i.e. when tand = Evi/v 2 - For the sake 

of brevity, we have not written down the interactions 
of the heavier Majorana mass eigenstates of the theory. 
Although such terms are only relevant when the mass 
splittings are not too large, we do include them in our 
final analysis. 

The analytic forms for the Higgs and gauge interac¬ 
tions we have written down so far are not completely 
general, since they depend on the sign of the mass terms 
for the physical states of the theory. In particular, if 
upon diagonalizing to the mass eigenstate basis we en¬ 
counter states with a “wrong sign” mass term, a field 
redefinition that preserves the Majorana character of the 
field must be performed, e.g. x This transfor¬ 

mation preserves the canonical form of the kinetic terms, 
but may introduce additional y^’s or Vs in the interaction 
terms. These changes lead to physical consequences. The 
explicit analytic expressions written down in this paper 
must then be appropriately modified if this is indeed the 
case. 


B. General Constraints 

To be considered as a realistic DM candidate, x must 
satisfy a handful of experimental constraints: 

1. The thermal relic density of x must satisfy 

« 0.1199 ±0.0027, in agreement with current 
measurements from WMAP and Planck [1]. Fur¬ 
thermore, the annihilation rate must lie below the 
95 % CL upper limits from gamma-ray observations 
of dwarf spheroidal satellite galaxies of the Milky 
Way [67]. 

2. The scattering rate of x with nuclei is below the 
current spin-independent and spin-dependent lim¬ 
its from LUX [63, 68]. 

3. If rux < m.h/2, then BF(h —>■ xx) 7 19%, com¬ 
ing from global fits to the observed Higgs couplings 
[69]. Note that our results do not change appre¬ 
ciably if we impose weaker bounds on the invisible 
Higgs decay at the level of 50 — 60%, as obtained 
from direct searches [54, 70]. 

4. If rUx < mzl‘2,, then V{Z -A xx) 7 2 MeV, coming 
from LEP precision electroweak measurements near 
the Z-pole [65]. 
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model=du, Ms=l50 GeV, Md=500 GeV 


model=dd, Ms=150 GeV, Md=500 GeV 




FIG. 3. Spin-independent tree-level direct detection cross section normalized to the LUX excluded cross section as a function 
of tan0, for several values of y and mu, in the {left) du and {right) dd singlet-doublet models. The default parameters are 
Ms = 150 GeV, Md = 500 GeV, and tan (3 = 5. 


5. No other constraints from the LHC or other col¬ 
liders are violated. In particular, we will consider 
LEP and LHC direct searches for a heavy Higgs 
and electroweakinos, and mono-6 constraints. 


We will represent the regions with the correct relic 
density in black in our summary plots of Figs. 4 and 5. 
Regions where the annihilation rate is above the 95 % 
CL upper limits from gamma-ray observations of dwarf 
spheroidals will be represented in purple. 


Relic Density. For DM masses < 80 GeV, the 
dominant annihilation channels governing freeze-out and 
annihilation today are s-channel exchange of Z and/or 
A (if A is sufficiently light). Since x is Majorana, for Z- 
exchange the s-wave contribution is chirality suppressed 
by the mass of the final state SM fermions. Meanwhile, 
H-exchange is also suppressed by fermion mass but can 
be enhanced for large tan ,5. 

For more massive DM, specifically when y is suffi¬ 
ciently heavier than mw, x can annihilate into pairs of 
W±’s or Z’s, and both of these processes are in gen¬ 
eral unsuppressed by the y relative velocity. Once 
is taken to be larger than a few hundred GeV, annihila¬ 
tions to pairs of scalars become relevant, which depends 
on additional couplings in the full Higgs sector. We will 
therefore restrict our study to the parameter space where 
y is sufficiently light such that final states including iL, 
A, and do not contribute significantly to the calcu¬ 
lation of . With this simplifying assumption, we can 
safely ignore the heavy Higgs self-interactions present in 
the full 2HDM Lagrangian. 

We find that in the calculation of resonances 

and CO- annihilations can play an important role in set¬ 
ting the relic abundance in certain regions of parameter 
space. Since a proper calculation of flh^ requires a care¬ 
ful treatment of such effects, we implement our model 
with FeynRules 2.0 [71] and micrOMEGAS_3.6.9 [72] in 
the calculation of the relic abundance and numerically 
scan over the parameter space of our different models. 
In doing so, we have checked the output of the dominant 
annihilation channels analytically. 


Direct Detection. The dominant contribution to the 
SI scattering rate of y off of nuclei is from exchange of the 
CP-even scalars h, H. The SI y scattering cross section 
per nucleon is 


SI, per nucleon 
^0 


-Z^/prot + {A — Z)fneut' ^ 
TT L A 


, (23) 


where A and Z are the atomic mass and atomic number, 
respectively, of the target nucleus, and is the re¬ 
duced mass of the WIMP-nucleon system. /„ (where the 
nucleon n is a proton or neutron) is the effective WIMP 
coupling to nucleons, which can be written in terms of 
the quark couplings as 


fn — 


E 3 i _ f (") I _ f (") 

/t, + 27^™ ^ 

■ q—u,d,s ^ 

where aqjrrig is defined to be 


q—c,b,t ^ 


ujq 

rrin 


q0H = 


2 ' 2 
m/ mjj 

cot /3, if g = up-type 
— tan 13, if g = down-type . 


(24) 


(25) 


In Eq. (24), we take the quark mass fractions to have 
the values = 0.02, = 0.026, = 0.014, 

= 0.036, = 0.043. By definition, 

= 1— The effective WIMP-Higgs 

q—u,d,s 

Yukawa couplings in Eq. (25) are the same ones given 



















explicitly in Eqs. (B6)-(B9). Throughout, we demand 
that the rate given in Eq. (23) is below the SI constraints 
from LUX [63]. We will represent this constraint in red 
in our summary plots of Eigs. 4 and 5. 

As we will see later in Sec. VC, the sign of tanO can 
have important effects on the direct detection rate. Al¬ 
though it will only slightly alter the DM annihilation rate 
throughout its thermal history, a negative tan 9 can al¬ 
low for a suppressed SI scattering rate as seen by LUX. 
To illustrate this effect, we show in Eig. 3 the SI nucleon 
scattering rate (normalized by the LUX limit) for the 
du and dd models as a function of tan0 and for various 
representative choices of y, Ms, and Md- 

In the left plot of Fig. 3, we show the normalized rate 
for the du model. For tanO < 0, Eq. (B6) and Eq. (25) 
imply that the effective light and heavy Higgs couplings 
with nucleons can be comparable and opposite in sign. 
As can be seen in Fig. 3, this destructive interference be¬ 
tween h and H exchange is a generic feature of the du 
model for tan 9 < 0, however the exact point of max¬ 
imum cancellation depends on the specific choice of y, 
Ms, Mo, and tan/?. Furthermore, independent of the 
chosen benchmark of couplings, for large tan/?, there is 
always a suppression at tan0 = yily\ =0, which can be 
understood as the point of suppressed mixing since both 
mixing terms yiVi = yVdCOs9 and y 2 V 2 = yvu sin 9 are 
small here. 

In the right plot of Fig. 3, we show the normalized SI 
scattering rate for the dd model. In this case, the rate 
can also be suppressed for tan 9 < 0. This is because 
both couplings X^h and X^h vanish at a blind spot of 
the theory when rriy, + Mo sin 29 = 0. As can be seen in 
the figure, the position of the blind spot does not depend 
much on the values of y and mo- For the Ms and Mo 
parameters in Fig. 3, this occurs at tan9 « —0.15, —6.51. 

We next consider elastic spin-dependent (SD) scatter¬ 
ing of X with nuclei via Z exchange, requiring that the 
DM-neutron cross section® is consistent with the recast 
constraints from LUX [68]. The SD scattering of x 
neutron is 


SD, per neutron 
^0 


12 ^ 


2 

X,neut 


( H “9^9 


neut) 


q—u,d,s 


, (26) 


where aq is the effective coupling with quarks, 
„ _ 9xZ9qa 


9qa = ) , if g = up/down-type . (27) 


In Eq. (26), we take the quark spin fractions 
to be = -0.42, = 0.85, and 

^(neut) _ _Q Qg jg coupling of DM with 


® The constraints on SD scattering with neutrons is generally 
stronger than that coming from protons. 


the Z boson and is given in Eq. (22). Here, e is the elec¬ 
tric charge of the electron (e > 0), and tw is the tangent 
of the Weinberg angle. 

The constraints on the SD scattering off of neutrons 
from LUX become more important as the doublet mass 
Mo is decreased and the doublet fraction is correspond¬ 
ingly enhanced (of course, if y is purely doublet, then 
g^z vanishes completely as seen in Eq. (22)). We will 
represent this constraint in orange in our summary plots 
of Figs. 4 and 5. 

Invisible Decays. The constraints from the invisi¬ 
ble widths of the SM Higgs h and Z are relevant when¬ 
ever < mh/2, 771^/2, respectively. The invisible 
branching fraction of the Higgs is constrained to satisfy 
BF{h —>■ xx) ^ 0-19 at 95% CL, which comes from a 
global fit to the visible Higgs channels in which the visible 
Higgs couplings are fixed to their SM values [69] . Assum¬ 
ing that r{h —>■ SM) = r®^(/i —)■ SM), it follows that the 
constraint on the invisible width of the Higgs is approxi¬ 
mately r{h xx) ^ 0-99 MeV. From the Lagrangian of 
Eq. (19), the Higgs width into a pair of Majorana y’s is 
found to be 


r(h^ xx) = -^ruh 



(28) 


Similarly, electroweak precision measurements at LEP 
constrain r(Z —>• xx) ^ 2 MeV [65]. For < rnzl‘^, 
the Z width into a pair of x’s is 


r(^ ^ xx) = (1 



(29) 


We will represent the constraint from the Higgs {Z) in¬ 
visible width in gray (brown) in our summary plot of 
Fig. 5. 


Direct searches for new particles. The model con¬ 
tains additional (possibly light) new particles that can be 
looked for directly at the LHC: two neutral {H,A) and 
one charged Higgs boson {H^) and two neutral (x 2 , 3 ) and 
one charged {E) fermion, in addition to the DM candi¬ 
date, X = Xi- In Sec. HI A, we have already discussed 
the LHC bounds on neutral heavy Higgs bosons. Here 
we comment on the constraints on charged Higgs bosons, 
as well as on new fermions. 

The most relevant LHC charged Higgs searches are 
for the process pp —>■ , tbH^, with subsequent 

—)• Tiy [76, 77] or —)• tb [78]. These searches, 

performed with 8 TeV data, probe charged Higgs bosons 
in the multi-hundred GeV range, but only for very large 
values of tan/? (> 30 — 40), for which the coupling 
entering the charged Higgs production is en¬ 
hanced. In addition to direct searches, flavor transitions 
such as & —>■ S 7 can set interesting (indirect) bounds on 
the charged Higgs mass: in a Type-H 2HDM, charged 
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model = du. y = 0.5, tan0 = -1, tan^ = 7, = 300 GeV model = dd. y = 1.5, tan^ = -10. tar^S = 3, = ttih = ttih ^ = 300 GeV 



FIG. 4. Summary plot for the du model for y — 0.5, tan0 = —1, and tan/3 = 7 {left) and for the dd model for y = 1.5, 
tanS = —10, and tan/3 = 3 {right). Shown are contours of ~ 0.12 (black) and (blue-dashed) and exclusion regions 
from direct detection searches at LUX for the spin-independent process (red) and spin-dependent process (orange). We also 
show regions that are above the upper limits from gamma-ray observations of dwarf galaxies (purple). In both cases, the masses 
of the heavy scalars are fixed to be 300 GeV. 


Higgs bosons cannot be lighter than around (300-400) 
GeV [79] . In the following, we will always fix the charged 
Higgs mass to 300 GeV and tan/3 < 10 in such a way as 
to avoid constraints from flavor transitions and direct 
collider searches for heavy Higgs bosons. 

LHC searches for electroweak Drell-Yan production 
can set interesting bounds for the new fermions arising in 
our model. In particular, the model contains one charged 
and two neutral fermions, in addition to the DM. These 
fermions are produced either in pairs through a Z bo¬ 
son {pp —>■ XiXi) where Xi = X2,3j E), or in associated 
neutral-charged production, with the exchange of a IT 
boson {pp —>• X 2 , 3 E). Generically, the latter produc¬ 
tion mode has the most relevant LHC constraints. LHC 
searches for supersymmetric Wino associated production 
giving a 3£ -b (or 2£ -b) MET signature can already set 
bounds on the Wino mass at around 400 GeV for mass¬ 
less lightest supersymmetric particles (LSPs) [80-82]. In 
Sec. VD, we will discuss how this bound can be inter¬ 
preted in our model. 


C. Discussion on DM above 100 GeV 

Among all of the interactions of our DM candidate 
X, the coupling to the pseudoscalar A is key in opening 
up viable parameter space. Annihilation through an s- 
channel pseudoscalar is s-wave, and if it happens not too 


far from the pseudoscalar resonance, 

l-4m2/m^<0.5, (30) 

then it is possible to obtain ~ 0.12 even for a large 
singlet fraction. And unlike the interactions of y with 
h and H, y scattering off of quarks via exchange of the 
pseudoscalar A is spin-dependent (SD) and further kine¬ 
matically suppressed by four powers of the momentum 
transfer, g'*, where q is typically of order 100 MeV. 

In exploring the parameter space of the model, we have 
found that it is difficult to find thermal DM candidates 
very far away from pseudoscalar resonances. When the 
relic density is not governed by resonances (or coanni¬ 
hilations) to any significant degree, there must be large 
mixing in the gauge eigenstate makeup of y in order to 
obtain the correct abundance. Since this mixing hinges 
on the Yukawa interactions (see Eq. (18)), this will also 
increase the DM scattering rate off of quarks. Such re¬ 
gions are generally ruled out by LUX, except near special 
blind spots, as discussed above. 

We present benchmarks for two of the models of Ta¬ 
ble H in Fig. 4. As previously mentioned, the du and ud 
models can be related to each other by replacing tan0 
with its inverse, and the uu model in the large tan /3 
limit reduces to only coupling the SM Higgs to the dark 
sector. We therefore only explore the parameter space 
for the models du and dd. 

In this section, we focus on > 100 GeV, which also 
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easily avoids constraints from the invisible width of the Z 
or h. Prospects for lighter DM {rriy^ < 100 GeV) will be 
presented below in Sec. VD. We also choose to work with 
sufficiently light y, such that annihilations to final states 
including one or more heavy Higgs are negligible. This 
therefore favors scalars of mass around a few hundred 
GeV in order for a pseudoscalar resonance to be relevant, 
and so we fix itia = mu = = 300 GeV in order to 

simplify the scan of the parameter space. We note the 
model would work just as well for heavier scalars and 
correspondingly heavier DM. 

In. Fig. 4, we show constraints in the (Md, Ms) plane 
for representative choices of y, tan 0, and tan /3 that 
give viable parameter space for thermal relic DM. Our 
choice of parameters satisfies the constraints from direct 
searches for heavy Higgs particles in the brT and bbb fi¬ 
nal states, as discussed in Section HI. Furthermore, since 
the pseudoscalar-DM coupling < 0(0.1) in the full 
Mjj — Ms plane we present, mono-5 searches are much 
less sensitive to this model. 

In both du and dd models, we clearly see the least con¬ 
strained region for thermal relic DM is where < 150 
GeV (m^ at or slightly above 150 GeV is in slight tension 
with gamma-ray observations of dwarf galaxies). The 
thermal relic line extends to larger Mo (larger singlet 
fraction) as the DM mass approaches the resonant region, 
which is centered slightly below = m^/2. This can 
be understood as thermal broadening of the resonance 
near freeze-out. Another feature in the thermal relic line 
can emerge if tan/3 is not too large: in the right frame of 
Fig. 4, for > 175 GeV, there can be dominant anni¬ 
hilations to top quarks through s-channel exchange of a 
light or heavy Higgs. 

Although the pseudoscalar resonance region will re¬ 
main viable for many different parameters, the proxim¬ 
ity or tuning of to mA/‘2 depends on the choice of 
tand or y. The thermal relic line in the figures will shift 
to larger values of Mo for a fixed Ms for large |tan0|. 
(The relic abundance is only slightly affected by the sign 
of tanO.) This is because, in both du and dd models, for 
large |tan0|, \y 2 V 2 \ |yi'(^i|: which then implies (using 

Eq. (18)) |A^f| Hence, g^z is largely unsup¬ 

pressed, as seen by Eq. (22), and to compensate, the 
overall singlet fraction must be increased by slightly 
decoupling Mo- Similarly, the doublet fraction of x is 
proportional to y, and so increasing y in either of the 
scans of Fig. 4 will generically shift the thermal relic line 
to larger values of Mo for a given value of Ms- 

Direct detection constraints are most relevant at lower 
Mo, where the DM singlet fraction is lower (see the red 
region of Fig. 4). Fig. 4 also illustrates various blind- 
spots in the SI direct detection rate. In the du model, for 
Mo ~ 100 — 200 GeV, the choice of tan0 < 0 suppresses 
nucleon scattering. In particular, from Eq. (B6), in the 
large tan /? limit the relative strength of the two different 
GP-even Higgs couplings is ^ih/Kn ~ (m^/Mr,) tan 6». 
Then, from Eq. (25), when tan0 = —1 the couplings 
of nucleons to h and H partially cancel, explaining the 


feature in the drop off in scattering rate for Mo ~ 200 
GeV and Ms ^ 100 GeV. In the dd model, shown in 
the right frame of Fig. 4, x-nucleon SI scattering is also 
near a blind-spot of the model where both X^h, X^h are 
suppressed. Here, fixing tan0 = —10, X^h ~ X^h ~ 0 
approximately when Mo/m^ « 5. 

The SD constraints (the orange region in Fig. 4) do 
not rule out much of the viable thermal relic parameter 
space that SI constraints do not already exclude. The 
importance of considering SD scattering via Z exchange 
is still illustrated in the du model, where spin-dependent 
limits are more powerful than spin-independent limits for 
Mo values of < 200 GeV. The SD limits do not depend 
on Higgs couplings and are able to constrain the parts 
of the parameter space close to where h and H exchange 
interfere and suppress the SI scattering rate. 

Finally, as already introduced in the previous sub¬ 
section, additional constraints might come from the 
LHG direct search of electroweak Drell-Yan produc¬ 
tion. However, in the regime with not too light DM, 
^ (100 — 150) GeV, there are no bounds even for 
Mo of around 200 GeV [80, 81]. 


D. Light DM and the GCE 

In this section, we investigate the viability of the mod¬ 
els to describe DM with mass below ^ 100 GeV. Al¬ 
though much of the physics near the pseudoscalar reso¬ 
nance is similar to that of the previous section, we addi¬ 
tionally require that the model could provide a reason¬ 
able fit of the Galactic Genter Excess. For model building 
in a similar direction that can additionally describe the 
3.55 keV X-ray line, see [83]. 

One simplified model that has received much atten¬ 
tion in its ability to describe this signal is just that 
of Eq. (1). Since the annihilation is s-wave, the rate 
can still be large today. Moreover, a relatively light 
pseudoscalar is favored in UV-complete realizations to 
get around numerous constraints [32-34]. The same is 
true in our models, since a pseudoscalar of around 100- 
200 GeV will be needed for the DM to annihilate near¬ 
resonance. As discussed in detail in Sec. H, we therefore 
implement the freedom in a general 2HDM to have a 
sizable splitting between the pseudoscalar mass and the 
heavy/charged Higgs mass and fix ruA = 160 GeV and 
■niH = = 300 GeV. This allows x to be relatively 

close to a pseudoscalar resonance, while the other scalars 
are heavy to evade direct detection and other constraints. 

We now add the following criteria to the list of de¬ 
mands enumerated in Sec. VB for x to be considered a 
realistic DM candidate consistent with the GCE: 

1. Since annihilations through a light pseudoscalar 
proceed dominantly to final state bottom quarks, 
we restrict the mass of x to li® the range 
50 GeV < < 70 GeV in order to fit the 
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FIG. 5. Examples of viable parameter space for the GCE. In the left panel we show the du model; in the right the dd model. A 
DM with the correct relic density, mass, and s-wave annihilation cross section for the GGE is found when there is an overlap of 
the solid black curve, the blue dashed region, and the green dashed region. Regions are excluded by direct detection searches 
at LUX for the spin-independent process (red) and spin-dependent process (orange). In brown and gray we show the regions 
excluded by Z —>■ invisible and by h —^-invisible, respectively (these place no constraints on the dd model parameter space 
shown). We also show regions that are above the upper limits from gamma-ray observations of dwarf galaxies (purple). In 
both models, the masses of the heavy scalar and charged Higgs are fixed to be 300 GeV. The pseudoscalar mass is 160 GeV. 


spectral shape of the GCE spectrum (this is repre¬ 
sented in blue in Fig. 5).® 

2. In order to ensure the approximate normalization 
for the GCE signal, the annihilation rate must 
satisfy 0.5 x 10“^® cm^/s < limt,_j,o (ctu) < 4 x 
10“^® cm^/s [30] (this is represented in green in 
Fig. 5). 

Limits on DM annihilation from AMS-02 observations of 
the anti-proton fraction [84, 85] are very similar to those 
from the dwarf galaxies and have additional astrophysical 
uncertainties, so for simplicity we do not include them. 

In Fig. 5 we show benchmark scenarios for ^100 
GeV, highlighting the regions that ht the GCE. We 
again consider the du and dd models and scan over the 
{Md,Ms) plane. We make a similar choice for y, tan9, 
and tan/3 as in Sec. V C. For both models, the qualitative 
behavior of the direct detection constraints is similar to 


The exact mass range that is preferred is dependent on system- 
atics. For annihilations to bb, the spectral shape of the observed 
emission has been found to be well fit by DM of mass as low as 
30 GeV and as high as 70 GeV [20, 21, 29, 30]. As a benchmark, 
we choose the upper half of this range since it is more viable for 
the model at hand. 


the previous section. Again, since the pseudoscalar-DM 
coupling ^ G>(0.1), mono-& searches would have very 
little sensitivity to the relevant parameter space. 

For both models in Fig. 5, the GCE excess can 
be explained while avoiding direct detection limits if 
X is mostly singlet-like and near a pseudoscalar res¬ 
onance. Constraints from gamma-ray observations of 
dwarf spheroidals may also be avoided if y is slightly be¬ 
low the pseudoscalar resonance, due to thermal averaging 
of the annihilation cross section. We also see the effects 
of annihilation near the Z, h poles (the latter only for the 
du model), which is visible in the thermal relic contour 
at around Ms ~ 45,60 GeV, respectively. This effect is 
not present for lim„_j.o (cru); for Majorana fermions anni¬ 
hilating through an s-channel vector mediator (or Dirac 
fermions with only axial couplings) there is no resonant 
enhancement in the s-wave contribution to av [86] and 
annihilation through an s-channel scalar mediator is p- 
wave for fermionic DM. 

The favored region for the du model is for 
Md ~ 300 GeV - 1 TeV, and - Mg = 60 - 70 GeV. 
Due to the relatively smaller mixing induced in the dd 
model, the appropriate GCE parameter space requires 
a smaller Mjj (200 — 500) GeV. We also emphasize 
that, although we have presented only two benchmark 
scenarios here, these choices of models and parameters of 
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the DM sector are not particularly special or highly fine- 
tuned. In particular, all models with a pseudoscalar that 
is not too heavy, tan 0 < 0, and not too small tan /3 could 
give a good fit of the GCE. Of course, pseudoscalars even 
lighter than 160 GeV would be suitable to obtain a large 
enough annihilation rate. For example, a pseudoscalar 
mass around 100 GeV is even better suited for the GGE. 
For lighter uia ^ 100 GeV, then there are currently no 
LHC direct searches for pseudoscalars (although this may 
be possible in the future [87]) and only very weak con¬ 
straints from heavy Higgs searches at LEP [88] are appli¬ 
cable. However, such very light pseudoscalars are more 
difficult to be achieved within our 2HDM scalar sector 
(see Appendix A for more details). 

We finally comment on the additional constraints we 
have to consider for these models with light DM, and con¬ 
sequently with new relatively light electrically charged 
degrees of freedom {E in Eq. (19) and Eq. (21)). We 
first note that the constraint from the invisible width of 
the SM h and Z are not strong and can exclude only a 
small region of parameter space in the du model at light 
Md, which is already excluded by LUX constraints. 

Additionally, having new light and SU{2) charged de¬ 
grees of freedom can introduce corrections to electroweak 
precision observables. However, we have checked that 
the contributions to the T parameter from loops of 
the new fermions is negligible (at most at the level of 
AT - 10-4 - 10-3 in the region of parameter space fa¬ 
vored by the GCE) [89]. 

Finally, constraints on the parameter space come from 
LHC direct searches for Drell-Yan production of elec¬ 
troweak particles. For the models shown in Fig 5, the 
electroweak spectrum needed to fit the GCE contains sev¬ 
eral new light fermions in addition to DM. In particular, 
for the dd model, we have two additional neutral fermions 
X 2,3 and one charged fermion E, all with mass close to 
Md ~ (200 — 500) GeV range and with splittings smaller 
than a few GeV amongst the states. Constraints from 
LHC searches of neutral-charged Wino associated pro¬ 
duction, resulting in a 3^ -I- (or 2^ -I-) MET signature, are 
the most important to constrain our scenario. In particu¬ 
lar, combining the 21 and 3£ searches, the ATLAS collab¬ 
oration sets a bound at uiwino 400 GeV, under the as¬ 
sumption of 100% branching ratio W^Z+^t 

and massless LSP [81]. This corresponds to an exclusion 
on a{pp —)■ —>■ W^Z -|- .^t) ^ 30 fb. In the 

following, we discuss how to interpret this constraint in 
terms of our model. 

As discussed in Sec. V, the fermion content of our 
model resembles the one of the MSSM with a Bino-like 
LSP and Higgsino-like NLSPs. However, in our model 
the heaviest fermions xi Xs have a sizable branching 
ratio X 2,3 —>■ Axi, as long as it is kinematically accessible 
{Mjj > 220 GeV for the benchmarks in Fig. 5). Then it is 
easy to check that, in the dd model, the cross section for 
PP X 2 , 3 E —>■ WZxiXi is always smaller than the ex¬ 
cluded cross section (30 fb) in the entire region of param¬ 
eter space for Mjj > 220 GeV. Below 220 GeV, the decay 


into a pseudoscalar is not accessible and the branching 
ratio for the decay X 2,3 —>■ hxi is not large enough to sup¬ 
press sufficiently the X 2 , 3 E —?> VFZxiXi channel. There¬ 
fore, in the dd model, the region Mu < 220 GeV of Fig. 5 
has already been probed by the LHC direct searches of 
Drell-Yan production of electroweak particles. On the 
other hand, the region of parameter space favored by the 
GCE in the du model has not been probed by these LHC 
searches yet, since the additional fermions are heavier 
(Md > 300 GeV). 


E. Future tests of the model. 

Upcoming spin-independent direct detection results 
can probe much of the parameter space for our model. 
For the du model, which has a larger SI rate compared to 
the dd model, future LUX data can test the entire GC ex¬ 
cess region and much of the region for heavy dark matter. 
The next generation of ton-scale direct detection experi¬ 
ments, such as LZ [90] and XENONIT [91], will be vastly 
more constraining. The LZ experiment can cover the en¬ 
tirety of the parameter space shown for the du model in 
Figs. 4-5, while for the dd model a significant portion of 
the parameter space can be reached, including the entire 
GC excess region shown in Fig. 5. 

Our model can be further probed at LHC Run H by the 
search of the various light degrees of freedom. In partic¬ 
ular, as shown in Fig. 2, pseudoscalar searches (with the 
pseudoscalar decaying either to taus or invisibly) will be 
able to test almost entirely the region of parameter space 
able to predict a thermal DM candidate, if both the DM 
and pseudoscalar are relatively light ( ttidm < /2 and 

niA below ~ 200 GeV). Furthermore, at the Run H of the 
LHC, with 300 fb-^ data, searches for Drell-Yan produc¬ 
tion of electroweak particles will be able to probe Wino 
masses as high as ^ 840 GeV [92], under the assumption 
of 100% decay for -A W +^t and for -A Z+^t- 
This can be translated into a bound on Md at the level 
of ~ 300 GeV in the dd model for the GCE region^. 

Finally, our model also predicts additional signatures 
that can be searched for at the LHC. As discussed 
above, Drell-Yan production can dominantly lead to 
PP X2,3E -a VFAxiXi) with A -A bb,TT, or even 
XiXi- Similar signals have A replaced by h. In par¬ 
ticular, for Md < 300 GeV, the new scalars will have 
sizable branching ratios into the Xi and E fermions: the 
charged Higgs can decay to DM and to a charged fermion 
A, E[^ -A Exi, producing a W +^t signature; the heavy 
scalar can decay to DM and to an additional neutral 
fermion H -A XiX 2 , 3 ) with X 2,3 possibly decaying to xi 
plus a Z, h, or A. It will be interesting to investigate 


^ Note that for Md above ~ 350 GeV, a new decay mode for E 
becomes available E —t suppressing even more the cross 

section for the process pp -A X 2 , 3 E -A WZ + 
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the potential of the pp —)■ tbH ^, —)■ W + and 

pp —)■ tt {bb) H, H ^ XiX 2,3 2^ (2&, 2t) + channels 
in probing our model. 

VI. CONCLUSIONS 

Viable theories of weak scale dark matter are under 
increasing pressure from the unyielding progress of cur¬ 
rent direct detection and collider experiments. Simplified 
models of dark matter that introduce preferred pseu¬ 
doscalar interactions with quarks explain the null re¬ 
sults of LUX, while still motivating interesting collider 
searches in the form of mono -6 or tt + processes. If 
the pseudoscalar is introduced in the context of a general 
2HDM, we have shown how the freedom in the scalar po¬ 
tential allow the pseudoscalar to be slightly decoupled, 
rriA < rriH Motivated by this framework, we 

have presented current and projected limits from missing 
energy and heavy Higgs searches at the LHC at energies 
of 8 and 14 TeV. 

We have also introduced a concrete realization of a 
model that couples a Type-H 2HDM to a fermionic dark 
sector consisting of a singlet and a pair of SU{2) dou¬ 
blets. The lightest Z 2 odd particle, y, is a dark mat¬ 
ter candidate that possesses couplings to the scalars af¬ 
ter mixing between the singlet and doublet. This model 
has analogues in the MSSM with bino-higgsino DM, and 
singlino-higgsino DM in the NMSSM, but we allow ar¬ 
bitrary couplings and scalar spectra. To simplify the 
presentation, we consider a discrete subset of couplings 
of the SU{2) doublets to the two Higgs doublets. Con¬ 
sidering direct detection, h and Z invisible decay, and 
other collider constraints, we identify regions of param¬ 
eter space with thermal relic DM below a few hundred 
GeV. In all of the parameter space that we have pre¬ 
sented, thermal relic y remains unconstrained if its mass 
is in proximity to a pseudoscalar resonance. 

For < 100 GeV, this model can describe the excess 
gamma-rays coming from the GC as measured from the 


Fermi satellite. A good fit to the gamma-ray data is given 
by DM annihilation to 6 -quarks with DM mass in the 
range 50-70 GeV. For these masses, we find that the pre¬ 
ferred scalar spectra has a light pseudoscalar of around 
160 GeV and heavier scalars at around 300 GeV. Future 
direct detection data from LZ will be able to thoroughly 
test these scenarios. There are also numerous LHC tests 
of the model, in the context of the GC excess. Precision 
Higgs coupling measurements can constrain the 2HDM 
that give rise to such light pseudoscalars. LHC searches 
for jets plus MET or the (in)visible decays of new heavy 
scalars can cover much of the parameter space at the 
14 TeV LHC. Consistency with the GCE generically also 
predicts new electrically charged states in the mass range 
of 200 — 1000 GeV. At Run 2, direct production of these 
new states can be searched for in the M + I^t channel. 
In addition, there are additional new collider signals with 
MET, heavy flavor, and gauge or Higgs bosons that result 
from the model. These new signatures may be valuable 
to study in future LHG data, independently of the spe¬ 
cific model we have presented. 
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Appendix A: Higgs spectrum 

In this Appendix, we comment on a feasible configuration of quartic couplings that is able to produce the Higgs 
spectrum presented in Sec. VD. In particular, we want to achieve a sizable splitting between the pseudoscalar mass 
mA and the heavy neutral and charged Higgs, as well as an approximate alignment limit cos(a — /3) ~ 0, for which 
the 125 GeV Higgs (h) has SM-like properties. 

To achieve a large splitting between the pseudoscalar and the charged Higgs, we need a sizable A5 — A4 (see Eq. 
( 6 )). Additionally, the value of A 2 is fairly constrained by the requirement of having the lightest scalar with mass 
at around 125 GeV (see Eq. (9)). Furthermore, in the regime of sizable tan/3, to get the alignment limit, we need 
to have the off-diagonal term for the mass matrix in the {hu,hd) basis much smaller than the diagonal terms. This 
translates into the condition A 3 -I- A 4 ^ Finally, our quartic couplings have to satisfy the conditions for a 

potential bounded from below, and, in particular, A3 > — A4 -I- jAsj — ^/Ai^ [93]. Putting together these conditions, 
we learn that also the value of Ai should be sizable. 

A possible benchmark that produces the Higgs spectrum presented in Sec. VD is given by 

Ai = 2.2, A 2 = 0.24, A 3 = 1.4, A 4 = -1.05, A 5 = 1.05, tua = 160 GeV, (Al) 
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that leads to rrih ^ 126 GeV, tuh ~ w_f/± 300 GeV, and relatively close to the alignment limit (cos(a — /3) ^ 0.1). 

From this discussion, we also learn that it will be difficult to have a much larger splitting between the pseudoscalar 
and the charged/heavy Higgs, if we want to maintain a stable potential and not too huge values of the quartic 
couplings. 

This benchmark scenario should be compared to the quartic coupling configuration of the MSSM: 


^MSSM _ ^MSSM _ 9^ + 9l ^ Q ^MSSM _ 

AfssM =~- 0.21 , 


0.08 , 


(A2) 


where g and gi are the SU{2) and t/(l) gauge coupling, respectively. 


Appendix B: Higgs couplings in singlet-doublet model 

We will parametrize the Higgs doublets as 

$ =±( V2N+2H^ \ 

72 + ziVi^ 2 ^/ ’ 

where, in the alignment limit, the coefficients ^re 

Nj} = cos/3 , = sin^ , /Vf = sin/3 , /V^ = - cos/3 


(HI) 


= — sin /3 , /V^ = cos /3 , = — sin /3 , = cos /3 . 

The Higgs couplings of Eq. (19) in terms of mixing angles are 


(B2) 


Axft = {yiN^N^ + y2N}N^) 

Kh ^ -^iV| (yi/V^iVf + y^mN^) 

Ka ^ {y,N^N^ - y2N^N^) 

Ks = [yiNt + y23V2+) 

-y2N+) . (B3) 

Analytic results are obtained by taking partial derivatives with respect to Vi or V 2 of the characteristic equation 
Det(Mneutrai — ^x) = ^ and then solving for the couplings A^^i and A ^2 defined by 


_ 1 dm^ _ 1 dm^ 


(B4) 


With these unphysical couplings in hand, the physical ones to h, H, A can be obtained, 

Ax^i = X^i + N 2 A;^2 
A^/f = At^i + N2 A^2 

AxA = —Ni + N 2 X^2 , (B5) 


where these couplings are exactly the same ones as defined in Eq. (B3). By following this procedure for each of the 
dw, ud, dd, uu models, we find: 


A 


du _ 

Xh — 


1 2 (1 + COS 2/3 COS 20) + Md sin 2/3 sin 2d 

2 ^ 2M|, + 4Msm^ — 6m^ + \y‘^v‘^ (1 + cos2/3cos 20) 
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1 2 iTT-x 2/3 COS 20 — Mjy cos 2/3 sin 29 

2^ 2M‘^ + AMsm^ — 6m’^ + {I + cos2P cos29) 

1 2 mj(. sin2/3 + Md sin20 

2^ 2M|, + 4Msto^ — 6m^ + (1 + cos2/3cos26*) 


(B6) 


A 


ud 

Xh 


y^v 


(sin^ /3 cos^ 6 + cos^ /3 sin^ 0) + sin2/3 sin20 
2M'^ + AMsm^ — 6m‘^ + (1 — cos 2/3 cos 20) 


A 


ud 

XH 


1 2 sin 2 ^ COS 20 + Md COS 2/3 sin 20 

2 ^ 2M^ + AMsrriy^ — 6m‘^ + (1 — cos 2/3 cos 20 ) 


1 2 2/3 + Mu sin 20 

2^ 2M|, + AMsm^ — 6m^ + (1 — cos 2/3 cos 20) 


(B7) 


xdd _ 2 2 o _ + Mp sin 20 _ 

^ 2 M|, + iMsm^ — 6 m 2 + y 2 v 2 cos^ /3 

. dd _ 1 2 • QO + Mjisin20 _ 

2^ 2^1"^ + AMsm^ — 6rn^ + y'^v'^ cos^/3 

dd 1 2 • ofl cos 29 

xA 2^^® 2M|) + 4M5TO;^ — 6to^ + y^t>2 cos^/3 


(B8) 


\ UU _ 


2 -20 

y II sin p 


+ Mp sin 20 

2 M|, + AMsm^ — + y^n^ sin^ /3 


\ UU _ 

\h--2 


y^v sin 2/3 


+ Mp sin 20 

2M|) + AMsm^ — Gm^ + y2v2 sin^ /J 


•^xA = -;^0^wsin2/3 


mv cos zfy 


2Mjj + AMsm^ — Gm| + y2w2 sin^ /3 


(B9) 


The M(i model is related to the du model by0—>'7r/2 — 0 along with the pseudoscalar coupling, A;,^^, flipping sign. 
The UU model, instead, is related to the dd model by ,3 —>■ /3 + 7r/2. 
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